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A detailed analysis of the capacitance of gated GaAs/AlxGai_xAs heterostructures is presented. 
The nonlinear dependence of the capacitance on the gate voltage and in-plane magnetic field is 
discussed together with the capacitance quantum steps connected with a population of higher 2D 
gas subbands. The results of full self-consistent numerical calculations are compared to recent 
experimental data. 
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I ' In modulation doped GaAs/AlxGai_xAs heterostructures a charged channel is formed in GaAs near the interface 
, due to the electron transfer from donors in AlxGai_xAs. If a metal Schottky gate is evaporated on the AlxGai_xAs 
d • surface the electron density N of the GaAs inversion layer can be controlled by the gate voltage Vg . In strong inversion 
and at low temperatures the impurity depletion charges are essentially fixed and the charge induced by the variation 
of the gate voltage goes entirely into the inversion layer. In this regime the heterostructure behaves like a capacitor 
with one metal electrode and the second electrode represented by the channel in GaAs. 

The differential capacitance (per unit area) of a gated structure, C = d\e\N/dVg, is traditionally written in 

the form 

O 



O 

ON 



Here Cb = e/dh denotes the barrier capacitance determined by the AlxGai_xAs barrier width db and the dielectric 
constant e. The second term takes into account the distance z of the centroid of the inversion layer charge from the 
interface. The numerical prefactor 7 results from approximating (see e.g. ^ and Q) the dependence of the subband 
energy Eq on N by dE^/dN cx z. In GaAs/AlxGai_xAs, 7 lies in the interval 0.5 - 0.7. The last term of equation (|l|) 
' expresses the dependence of the capacitance on the density of states g of the 2D electron gas. 

Recently the capacitance of gated GaAs/AlxGai_xAs structures has been measured with a high accuracy by Hamp- 
^ , ton et al [||, both as a function of the gate voltage and of the in-plane magnetic field. The main results of the 
''O ' measurements are, with a kind permission of authors, reproduced in figure ^ The aim of this paper is to present 
, a detailed analysis of the data based on the full self-consistent numerical calculation, going beyond the validity of 
^ ■ equation (Q). 

' A band diagram of the charged gated structure is sketched in figure Boundary conditions for the electrostatic 
potential $(ia) = and d^/dz(La) = correspond to the assumption of fixed impurity depletion charges and to a 
charge neutrality of the whole structure. The gate voltage Vg is a difference between the chemical potential Ep/\e\ 
of the gate and the chemical potential Ep/lel of the inversion layer. At the metallic gate the chemical potential 
■ is determined essentially by the electrostatic potential ^g since the difference is a constant built-in voltage of the 
Schottky barrier Vb- We obtain from the Poisson equation 

^^^^\e\N + K. (2) 

The constant K is a contribution to the potential due to frozen out impurity charges and due to the discontinuity at 
the GaAs/AlxGai„xAs interface. Consequently, the gate voltage is given by 

Vg = ^ + ^!d:I\e\N + K + VB (3) 
|e| e 

and the reciprocal value of the capacitance reads 

= C^' + (4) 

where 

C^'-4^(-N + ^) ■ (5) 
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Thus, the gated structure capacitance can be again considered as a barrier capacitance Cb and a channel capacitance 
Cc in series. The reciprocal value of the channel capacitance is a sum of two terms: the electrostatic one, related to 
the concentration dependent position of the channel centroid in GaAs and the thermodynamical one originating from 
the concentration dependent chemical potential in the semiconductor. Both these quantities can also depend on the 
applied magnetic field. 

Up to now a very general analysis of the gated structure capacitance has been carried out regardless the detail 
properties of the inversion layer. In the next step we summarize basic formulae describing the layer electron structure. 

In a zero magnetic field the standard envelope function approximation yields the Schrodinger equation for a con- 
duction electron 

^^\e\Vr^onf{z)~E^ij{r)^0 (6) 

where the confining potential Vconf is a sum of the electrostatic potential $, obtained within the Hartree approx- 
imation, and the exchange-correlation term V^c, describing the effects of electron-electron interaction beyond the 
Hartree approximation. We use the sirnple analytic paramertization of Vxc within the local-density-functional formal- 
ism, suggested by Ruden and Dohler [^|. Thus, the self-consistent solution to coupled Schrodinger equation (||) and 
Poisson equation for the electrostatic potential (f> gives the electron structure of the conduction channel, including the 



many-body corrections to Eq and z [gj . 

Looking for the eigenfunction of the equation in the form 



V'(r) = exp{ikxX + ikyy)(j)n{z) (7) 



we obtain 2D subbands of allowed energies 



E{k,,ky,n)^—{k'^ + k') + E„. (8) 
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Equations (|^) and are the mathematical expression for the separability of the free in-plane {x, y-direction) and 
bound out-of-plane (z-direction) components of the electron motion. 

Firstly, let us assume only the lowest subband to be occupied. Then, due to the parabolic dispersion of the in-plane 
energy, the 2D density of states per unit area g is a universal constant m/Trh? and the 2D concentration of the electron 
layer reads 

N - g{EF - Eo) (9) 

where Eq is the lowest energy of the bound out-of-plane motion. 

The average position zq of the electron layer is determined by the z-dependent part of the wavefunction as 

zE^ J z\Mz)\^dz . (10) 

Now, the attention can be turned back to the heterostructure capacitance. Due to and ( p^ ) quasi-analytical 
relations for both the thermodynamical and electrostatic parts of the reciprocal channel capacitance are available. The 
electrostatic term is related exclusively to the out-of-plane component of the electron motion. In the thermodynamical 
term both components of the motion play a role. However, according to their influence is independent and we 
obtain the channel capacitance in the form 

I _z^ N dz^ ttTi^ 1 dEo 
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Different 2D gas concentrations (or gate voltages) directly yield different shapes of the confining potential Vconf 
and thus the bound state eigenenergy Eq is modified. For confining potentials which do not depart substantially 
from an exactly triangular well one obtains the scaling laws Iq cx N^^^^ and Eq oc TV^/'^ from them the relation 
dEQ/dN cx 'zq and the equation (|l|) is regained. The quantitative deviations can be expected for more realistic 
shape of the confining potential of standard gated structures and, in the case of gated symmetric wells or inverted 
structures, the above mentioned scaling law is not valid at all. Therefore both dEQ/dN and 1^ have to be determined 
independently in the course of the self-consistency procedure. 
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In magnetic fields parallel to the heterostructure interface the separability between the in-plane and out-of-plane 
components of the channel electron motion is lost. Assuming the magnetic field B || y the wavefunction (j){z) becomes 
kx dependent. As discussed in more details in this results in a charge redistribution of the 2D gas and in a 

modification of 2D subbands. The centroid zq of the channel first shifts to the bulk GaAs and then, after the magnetic 
field reaches a critical value, it returns back to the heterostructure interface. The deviation of the energy dispersion 
E{kx) from the parabolic shape yields a density of states which form depends both on the magnetic field and the 
electron concentration Finally, the subband edge is shifted to higher energies by the in-plane field. It follows 
from the above mentioned effects on the electron structure that the electrostatic term and both contributions to the 
thermodynamical term are modified by the in-plane magnetic field in a rather complicated way. 

A sufficient amount of information about the heterostructure capacitance has been given to open the last problem 
we want to deal with. Except the very general formula (||) previous conclusions are restricted to 2D electron systems 
with one occupied subband. Now we extend the discussion to multi-subband systems since the filling or depleting of 
the higher subband is reflected in the observed capacitance quantum steps. A heterostructure where the number of 
occupied subbands is changed by the gate voltage in the absence of the magnetic field will be considered. A discussion 
concerning the in-plane magnetic field induced depopulation of higher subbands would be analogous. 

In the case of two occupied subbands we obtain the channel capacitance 

I z N dz TiT? 1 dE 
2" " 7 ^ TrfTV ^ 2H^ ^ N^rf^ 

in the form similar to ( [III ) but with zq and i?o replaced by 

^ = ^ ' (13) 



E=^^. (14) 

The quantities A'o = (Ep — Eq) m/nh^ and A^i = {Ep — Ei) mj-Kf? are the concentrations of electrons in the first 
and second subband, Nq + Ni = N . 

Two terms in equations (pT|), ( p^ ) do not contain the derivatives with respect to N. The first one, proportional to 
the centroid coordinate, does not cause any step when (|ll] ) is replaced by (|lj) after reaching the critical concentration 
Nc since 'z{Nc + 0) = zo(A^c — 0). The step in the second term is a universal constant and reflect the fact that the 
density of states is doubled above N = Nc- The additional steps come from the terms with derivatives which can be 
expressed for N > Nc similarly as for the case of single subband occupancy with a help of dzl/dN , dEi/dN, i = 0,1. 

The previous analysis shows the principal reasons for the gate voltage or in-plane magnetic field dependency of the 
one-subband channel capacitance as well as for steps in capacitance when a higher subband is populated. To obtain 
results comparable to experimental data full self-consistent numerical calculations, treating the electron-electron 
interaction beyond the Hartree approximation both in zero and non-zero magnetic fields, are unavoidable since the 
simple scaling rule leading to (|l|) is quite unrealistic in the presence of the magnetic field and for the double subband 
occupancy. 

Figure |l] presents data of low temperature capacitance measurements on an AlxGai_xAs/GaAs sample with x — 0.3 
and channel concentrations from about 1 x 10^^ to 4 x 10^^ cm~^. Since the shape of the confining potential depends 
also on depletion charges parameters which are not exactly known we performed numerical calculations for Na = 
4 X lO^^cm"'^ and several frozen- in depletions lengths La- Due to an unambiguous relation between the gate voltage 
dependency of the capacitance, including the position of the quantum step, and the shape of the quantum well it is 
possible to use the experimental curve measured at By = for determining La- This way we obtained La — 550 nm. 
(We do not discuss the height of the quantum step since a current leakage occurred at the gate voltage Vg = +150 mV, 
which might quantitatively modify the measured voltage dependency of the capacitance.) As shown in figure |^, if 
the theoretical and experimental gate voltage dependencies at the zero magnetic field are in a good agreement, the 
theory predicts, with a very reasonable accuracy, the behaviour of the capacitance for non zero in-plane magnetic 
fields. The impurity scattering of electrons together with the roughness of the interface are responsible for smearing 
of the quantum step, we simulate their influence by introducing the Dingle temperature To « IK. The effect of the 
magnetic field on both the electrostatic and thermodynamical terms in the reciprocal channel capacitance is shown 
in figure ^ to illustrate the violation of the formula (p. Note that the term 1/e^ dEp/dN cannot be simplified due to 
the field and concentration dependence of the density of states. 

In conclusion, we want to point out that the capacitance measurements provide a lot of information about 2D 
electron gas properties but the data have to be interpreted very carefully. It is well known that in structures shown 
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in figure ^ the 2D channel with higher concentration is confined closer to the GaAs/AlxGai_xAs interface resulting 
in a smaller distance between electrodes. However, this effect is not reflected directly in the capacitance due to 
simultaneous shift of bound states energies. Similarly, the charge redistribution by the in-plane magnetic field is 
accompanied with the shift of subband edges. Moreover, the modification in the density of states contributes to 
the magnetic field induced changes in the capacitance. Finally, the discontinuities in the 2D density of states are 
responsible only partially for the capacitance quantum steps. We have seen that also the ratio between the filling of 
subbands and the relative positions of electron layers in each subband play an important role. 
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FIG. 1. Experimental curves of the capacitance as a function of the gate voltage and of the in-plane magnetic field 
determined by Hampton et al Q. 

FIG. 2. Schematic band diagram for a gated GaAs/AlxGai_xAs heterostructure. 



FIG. 3. Self-consistently calculated gate voltage and in-plane magnetic field dependencies of the relative capacitance. For 
B = the results for three depletion lengths of acceptors are presented, Na — 4 x 10^* cm~"^. Magnetic field dependencies are 
plotted only for selected gate voltages and La = 550 nm to make comparison of experimental and theoretical curves easier. 



FIG. 4. An increment of the reciprocal channel capacitance 5C^ ^ — ^{B) — ^(0) (dotted line) in the magnetic field, 
decomposed in a thermodynamical term 1/e^ dSEp/dN and two contributions N/edSz/dN, 5'z/e to the electrostatic term. 
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